home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / sim2.f < prev    next >
Text File  |  2002-06-11  |  667b  |  29 lines
  1. C
  2. C PAGE 354-355: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: SIM2.FOR
  5. C
  6. C BUFFON'S NEEDLE PROBLEM SIMULATION (RANDOM)
  7. C
  8.       DATA  L,N,M/256,1000,0/ 
  9.       PI2 = 2.0*ATAN(1.0)   
  10.       ANS = 1.0/PI2 
  11.       DO 2 I = 1,N
  12.         W = RANDOM(L)       
  13.         V = PI2*RANDOM(L)   
  14.         IF(W .LE. SIN(V))  M = M + 1  
  15.         IF(MOD(I,100) .EQ. 0)  THEN   
  16.           PROB = REAL(M)/REAL(I)      
  17.           PRINT 3,I,PROB,ANS
  18.         END IF
  19.    2  CONTINUE    
  20.    3  FORMAT(5X,I5,2F10.5)
  21.       STOP
  22.       END 
  23.   
  24.       FUNCTION RANDOM(L)
  25.       L = MOD(16807*L,2147483647)
  26.       RANDOM = REAL(L)*4.6566128752458E-10      
  27.       RETURN
  28.       END 
  29.